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ABSTRACT 

The CP violating asymmetries for Cabibbo suppressed charged D meson decays in the 
standard model are estimated in the factorized approximation, using the two-loop effective 
hamiltonian and a model for final state interactions previously tested for Cabibbo allowed 
D decays. No new parameters are added. The predictions are larger than expected and 
not too far from the experimental possibilities. 
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It is well known that CP violating effects show up in a decay process only if the decay 
amplitude is the sum of two different parts, whose phases are made of a weak (Cabibbo 
- Kobayashi - Maskawa) and a strong (final state interaction) contribution. The weak 
contributions to the phases change sign when going to the CP-conjugate process, while the 
strong ones do not. Let us denote a generic decay amplitude of this type by 

A = A e iSl + B e iS2 (1) 

and the corresponding CP conjugate amplitude by 

A = A*e iSl + B*e l&2 . (2) 

The CP violating asymmetry in the decay rates will be therefore 

= |^l| 2 -|^| 2 = 2S(A£T) sinfo-fr) 

aCP |^| 2 + |^| 2 \A\ 2 + \B\ 2 + 2 $t(AB*) cos(5 2 -5 1 ) U 

Both factors in the numerator of eq. (|||) should be nonvanishing to have a nonzero effect. 
Moreover, to have a sizeable asymmetry the moduli of the two amplitudes A and B should 
not differ too much. 

While the phases of the weak amplitudes are expressed in terms of fundamental 
parameters of the theory (the angles and the CP violating phase in the CKM quark mixing 
matrix), the strong phases Si are in general unknown for heavy flavour decays. For this 
reason the asymmetries in charged B (or D) decays are less studied and less under control, 
theoretically @ , than the asymmetries in mixing assisted channels in neutral B decays . 
In a previous paper we have presented a rather successful phenomenological analysis of 
Cabibbo-allowed D decays. In that analysis a specific model for final state interactions 
was used, based on the assumption that the scattering phase-shifts are dominated by the 
nearest resonance. 

In this paper we exploit this model of final state interactions to estimate the CP 
violating asymmetries for two-body D decays in the standard model. 

These asymmetries have long be believed to be unmeasurably small. For Cabibbo 
allowed decays, this is a consequence of the smallness of the second interfering amplitude 
with different weak phase (B in ([!])), which can only be generated by mixing and doubly 
Cabibbo suppressed decay. In fact, asymmetries in D — > Kir have been recently suggested 
as a sensitive tool to observe possible "new physics" effects in CP violating amplitudes, 
taking advantage of the large strong phase differences. 

1 Although they are easier to measure experimentally, due to their self-tagging property. 
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In the Cabibbo suppressed case the penguin diagram contribution to the effective 
Hamiltonian, although its coefficient is rather small, provides a nonnegligible second am- 
plitude. We will therefore limit our calculations to Cabibbo suppressed decays. Moreover, 
only charged D will be considered. In fact, to describe final state interactions in D° decays, 
that receive contributions from not yet observed isoscalar resonances, new undetermined 
parameters would be needed. 

As we will see, the asymmetries are small but larger than previously expected 0, 
and in some cases they are 0(few 10 -3 ). To observe at 3 a level an asymmetry of 3 ■ 10 -3 
in a channel having a decay branching ratio of 0.5%, a number of charged D mesons equal 
to 2 ■ 10 8 is required, with 100% efficiencies. Present projects of tau-charm factories do not 
have a luminosity high enough to allow such a measurement, but it is not unconceivable 
that one could reach this sensitivity in a not too far future. 

The effective weak hamiltonian for Cabibbo suppressed nonleptonic decays of 
charmed particles is given by 

H ef f = ^§ V ud V: d [d Qf + C 2 Qi] + 

+ 5| v us v; s [d Ql + C 2 Q|] - (4) 



V2 

G p 

7f 



° F V ub V: b Y °i Qi + h.c. 



i=3 



In eq. (H) the operators are 
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7m (1 


-75)c/3 Y 
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(5) 



The operator Qf (Q|) i n ec L- (H) is obtained from Qf (Q 2 ) with the substitution 
(d — > s). In eq. (§) a and (3 are colour indices and in the "penguin" operators q (q) is to 



be summed over all active flavours (u, d, s). 
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We have evaluated the coefficients Cj using the two-loop anomalous dimension ma- 
trices recently calculated by A. Buras and collaborators [Q. The effect of next-to-leading 
corrections is particularly large for the coefficients of the penguin operators, as it has al- 
ready been noted in ref. in the case of AS = 1 decays. One problem that we have to 
face is the renormalization scheme dependence of the coefficients, that should be canceled 
by a corresponding dependence in the matrix elements of the operators. Unfortunately, 
we estimate the matrix elements of four-fermion operators using a specific model, and 
it is not evident to which scheme this corresponds @. We calculated the coefficients in 
both schemes (naive dimensional regularization and 't Hooft - Veltman scheme) for which 
anomalous dimension matrices have been given. We follow then the prescription of Buras 
et al., given in eq. (3.6) of ref. |§, to define "scheme independent" coefficients. Although 
not unique, this prescription is phenomenologically favoured, in that it reinforces the re- 
sults of the leading order calculation and provides large penguin contributions, that are 
welcome in the AS = 1 case. It turns out that this prescription gives for the dominant 
penguin operator coefficient Cq a value that is about two and a half times larger than the 
leading order prediction Cg . We stress that the results for Cq in both schemes are also 
considerably larger than C§ L . 

In ref. []2| we fitted the value of A% 1S , among other parameters, and obtained 
the result A^ s ~ 200 MeV. We therefore adopted this value in the calculation of the 
coefficients Cj. The results are reported in Table 1. 

Given the effective Hamiltonian, to obtain the decay amplitude 

A w (D^f) =<f\H eff \D> (6) 

we evaluate the matrix elements of the operators in the factorization approximation, ne- 
glecting moreover the colour-suppressed contributions, following || and 0. We also neglect 
the contribution of the dimension five operator 

Q= m c u<*cj^{l + l5 )ti p c?GA v , (7) 

that is generated at the two-loop level. Its Wilson coefficient has not been evaluated in 
||, where the masses of external quarks have been neglected. Its matrix element cannot 
be calculated in the factorized approximation and it would be anyhow suppressed by 1/N C 
in the large N c limit. 

2 An analogous problem for the dependence of the coefficients on the subtraction point, fi, 
already appears at leading order. 
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As an illustration, we discuss the explicit example D + — > K + K - a decay channel 
that turns out to have a sizeable asymmetry. The expressions for the matrix elements of 
the operators Qi (i=l,. . . 5) are easily obtained following ref. and are: 

< K+K*°\Q(\D+ > = < K+K*°\Q S 1 \D+ > = , 

< K+K*°\Qi\D+ > = -2M K * (e* • p K ) (m u + m d ) W PV , 

< K+K*°\Q S 2 \D+ > = 2 Mk* (e*- PK ) fx i _ ^/ M 2 > (8) 

< K + K*°\Q 4 \D+ > = < K+K* \Q*\D + > + < K+K*°\Q S 2 \D + > , 

< K+K*°\Q 3 \D+ >= < K+K*°\Q 5 \D+ >= . 

In eq. (||), in addition to self-explanatory symbols for particle masses and momenta and 
for the K* polarization vector e, two quantities need a definition: a cs ~ 0.8 is the axial 
charge and Wpy parameterizes the annihilation contribution. Wpy has been determined 
by a fit to the Cabibbo allowed decay rates in ref. @, with the result Wpy = 0.53 ± 0.13 
i. The operator Qq must first be written in a Fierz-rearranged form: 

Q 6 = -2 f(l - T o)cpu a (l + To )q a . (9) 

In the factorization approximation, its matrix elements are given in terms of matrix ele- 
ments of scalar and pseudoscalar densities. These in turn may be related to the matrix 
elements of the divergences of vector and axial currents, and thereby to form factors that 
have been discussed in H. The result is: 



< K+K*°\Q 6 \D+ > = AM K * (e*-p K ) 



M ° Wpy 



(m c + m d 

M 2 K f K a 



i (10) 



(to, + m u )(m c + m s ) 1 - M 2 K /M 2 Ds 

We note that the penguin operator Qq has, by far, the largest matrix element. 
Even if the coefficient Cg is considerably smaller than C2 or C\ its contribution would be 
dominant, were it not for the smallness of the CKM factor: \V ub V* b \/\V us V* s \ ~ 10" 3 . In 
fact, the Qq penguin operator gives the dominant contribution to the imaginary part of 
the amplitude A w , since all the CKM factors' imaginary parts are of the same order. 



Expressions analogous to (|8|) and (|10|) can be derived for the other two-body 



Cabibbo first forbidden decay channels of D + and 



3 We will assume the central values for the parameters determined in ref. [0 in the following 
estimates. 
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At this stage one would not have any CP violating asymmetry, since no strong 
phases from final state interactions have been included in the amplitudes A w (D — ► /). 
Rescattering, to be discussed in the next section, will provide the necessary nonzero strong 
phase-shifts. 



We make the assumption that the final state interactions are dominated by resonant 
contributions 0. 

In the mass region of pseudoscalar charmed particles there is evidence, albeit not 
very strong , for a J p = + resonance Kq (with mass 1950 MeV, width 201±86 
MeV and 52% branching ratio in Ktt 0), a J p = 0" #(1830) (with T = 250 MeV and 
an observed decay to K<j> @) and a J p = 0" tt(1770) with T = 310 MeV 0. These 



resonances should dominate, respectively, the rescattering effects for Cabibbo forbidden 
decays of — > PP, Df — ► PV and D + — ► PV, where P(V) indicates a pseudoscalar 
(vector) decay product. In || we assumed the existence of a J p = + resonance ao, that 
should be relevant for Cabibbo forbidden D + — * PP decays, with mass 1890 MeV and 
width ~ 200 MeV. We will not try to discuss the Cabibbo forbidden decay amplitudes for 
D° mesons; in this case several other unobserved isoscalar nonstrange resonances should 
determine the final state interactions and this would bring new unknown parameters in 
the calculations. 

The FSI effect modifies the amplitudes for Cabibbo suppressed D + and decays 
in the following way @: 

A(D -> V h P k ) = A W (D -> V h P k ) + c hk [exp(id 8 ) - 1] ^ c hlkl A W (D -> V^/ P fc/ ) = 

hlkl (11) 

= d^ fe A 27 + c hk exp(id 8 ) A 8 . 

In ( |TI| ) c^fc [d/ife] are the normalized c| fc = 1) couplings of a pseudoscalar (P) and a 
vector (V) octets to a pseudoscalar octet [27-plet]. The resonance P is present only in the 
octet case, and one has 

sin5 8 expUSs) = ^— ^ ^. (12) 

2{m~-m D )-iY{P) 

Similar expressions hold for A(D — > Ph Pk). 

The CP violating asymmetry of eq. (|3|) is given by: 

= 2 sm5 8 c hk d hk $S(A 27 A 8 *) 

a ° P c 2 hk \A 8 \ 2 + d 2 hk \A^\* + 2c hk d hk $l(A 27 A 8 *)cos5 8 [ J 
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To have a nonzero asymmetry the weak phases of A 8 and A 27 must be different. The 
phase difference is provided essentially by the penguin contributions, which dominate the 
imaginary parts of the amplitudes, and appear in A 8 only. On the other hand, it is true 
that in some channels there is no rescattering in our model and acp = 0. This is the case 
for instance for the decay amplitude A(D + — ► 7r + 7r°), which only contains the term A 27 . 

A particular treatment is needed for those decay channels where neutral kaons ap- 
pear. In fact, these particles will be mostly detected through their decays into two pions, 
and these decays are by themselves affected by CP violation (in the K system). One has 
therefore to be able to disentangle the CP violating effects in D and K decays. Define 
the "experimental" , time-integrated asymmetry for the D + decay to one charged and one 
neutral K meson: 

_ $q[T(D + -> K+K° -> K+n+n-)(t) - r(£T -> K~K° -> K-ir + ic-)(t)]dt 

Jq[T(D+ ->• K+K° -> K+n+n-)(t) + T(D- -> K~K° -> K-ir+ir-)(t)]dt 

(14) 

and recall the usual definitions of parameters for the neutral kaon decays 

AM =M Kl - M Ks , 

r=i(r s + r L ), (15) 

A(K L -> 7T + 7T") . 

V+-= ^ = e + e' . 

j4(AS — * 7T + 7T ) 

It is also convenient to define [ |i~4f : 

2 



/CO 



(am) 2 + r 2 



(r Kr7+_ + AM %r]+-) (1 - e" rT cos AM T) + 

+ (AMKfj + _ - r3jj + _)e" rT sinAMT 
Tl r 5 f(T) 



(16) 



^ (T) Tl (1 - e" r ^) + T S (1 - e" r ^) |t7 + _ | 2 ' 

7 ~ acp + 2 5ie . 
The resulting asymmetry is: 

A(T) = - 7 ~ ^> (17) 
1 - 7s(T) 

It is instructive to consider the limits of this expression for short (<< ts) and long (>> ts) 
time T: 

lim A(T) = acp 

(18) 

lim A(T) ~ a C p - 2 3fe . 

T — >oo 



Eq. (]T^) shows that to perform a measurement of the D + decay asymmetry in these 
channels will be essentially hopeless if the asymmetry acp turns out to be considerably 
less than 2 3? e, the asymmetry in Kl semileptonic decays. 

For D s decays, the equation analogous to (0) would be 



A (T) = J T [r(D+ -> TT + if -» Tr + TT + TT-Xt) - V {D~ -» „~ K° -> 7T~ TT + T^ ) (gjg 

J T [r(D+ -> 7T+K -> 7r+7r+7r-)(t) + r(D7 -> 7T-Z° -> n-n+n-)(t)}dt ' 



The resulting expression can be obtained from ([16]) and ([H]) by changing the signs of the 
terms containing e and rj^ Therefore 

lim A a (T) ~ a C p + 2 3?e. (20) 

T — >oo 



The decay amplitudes for Cabibbo-suppressed charm decays (and especially their 
CP violating asymmetries) depend on all the entries in the first two rows of the quark 



mixing matrix. In the Wolfenstein parametrization |T5j , extended to keep terms up to 
0(A 5 ) in the imaginary parts, the matrix is written as follows: 

/ 1-±A 2 A A\ 3 pe-* s \ 

V= -X(l + A 2 X 4 pe lS ) 1-|A 2 AX 2 (21) 

\AA 3 [l-(l-iA 2 )pe l<5 ] -AX 2 (1 + X 2 pe* s ) 1 / 

In (PH), A = 0.220 is the Cabibbo angle, we have chosen V c b = A X 2 = 0.046 and allowed 
different values for p, p = 0.50 ± 0.14, and for the CP violating phase 6. Following the 
analysis of |16|] and assuming for the top quark mass the "favourite" value (140 GeV), 



we considered both positive (0.47 -j- 0.81) and negative (—0.98 -. — 0.80) values for cos 5, 
corresponding to the two possible shapes of the unitarity triangle compatible with present 
data on B—B mixing and the value of the e parameter. The positive cos 8 values correspond 
to larger CP violating effects in K decays (e'/e and charged K decay asymmetries), in B 
decays and also in D decays. 

We have evaluated the branching ratios for decay into twelve different two-particle 
final states, and the corresponding CP violating asymmetries, both for D + and for Df. 
The results are reported in Tables 2 and 3, respectively. 

We note that the predicted branching ratios are in reasonable agreement with the 



existing data |1CJ], listed in the second columns, which have not been used in fitting the 
parameters. The variations in the theoretical predictions with p and cos 6 are less than 
10 -4 and therefore neglected in column three. 
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We also note that very recently a doubly Cabibbo-suppressed decay rate has been 
measured for the first time |T7[], with the result 



Mu^^i^.ir-. (22) 



With the parameters determined in our model predicts 

\ th . =3.2-10" 2 . (23) 



r(D+^<pK+)^ oo in _ 2 



T(D+ -> 07i , 

This ratio agrees within one standard deviation with the experimental result (^|) , however 
it should be noted that the predicted individual rates are smaller than the data by a 
somewhat larger amount, see Table 2. 

The errors in the asymmetries reported in columns four and five have been estimated 
by semi-dispersion in the numbers obtained varying one at a time p and cos 5 in the ranges 
previously mentioned. The asymmetries for an obtuse 5 (column five) are generally smaller 
by a factor two. There are some reasons to prefer the acute-angled solution for the unitarity 
triangle, though. Lattice QCD calculations ||18|| , as well as recent evaluations with QCD 



sum rules ||19|| , favour a rather large value for the B meson decay constant fs- This, 
together with the large top quark mass indicated by the electroweak data analyses pOfl , 
would select an acute-angled unitarity triangle. 

Hoping that nature is kind to us and cos 5 > 0, our results predict the largest 
asymmetries (~ 10 -2 ) for some D s rare decay channels. Considering however that D s are 
also difficult to produce, we consider the D + decay channels having the highest asymmetries 
(~ 3 • 10 -3 ) as more promising to observe CP violating effects. 

The best candidate is probably the decay D + — > K + K . Here the final state 
contains three charged particles 2/3 of the times, and the decay branching ratio has been 
already measured to be half a percent. These are the numbers we used in the discussion 
at the end of the introduction. A tau-charm factory at the ^"(3770) with an integrated 
luminosity C = 10 40 cm -2 would produce about 10 8 charged D mesons and therefore 
~ 1.5 • 10 5 decay events D+ -> K+K*° -> K + K~tx + and as many D~ -> K~K*° -> 
K~ K + n~ . Without including the further reduction due to actual experimental efficiencies, 
this number of events would imply for the asymmetry a statistical accuracy of 1.7 • 10 -3 . 
As we already said, the assumed luminosity is not enough to establish the effect, but at 
the same time not too far from this goal. 
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leading order 


naive dim. reg. 


't Hooft-Veltman 


"scheme indep." 


Ci 


-0.439 


-0.303 


-0.311 


-0.514 


c 2 


1.216 


1.139 


1.032 


1.270 


c 3 


0.0073 


0.0155 


0.0126 


0.0202 


c 4 


-0.0185 


-0.0398 


-0.0291 


-0.0460 




0.0057 


0.0099 


0.0084 


0.0127 


C 6 


-0.0210 


-0.0447 


-0.0320 


-0.0541 



TABLE 1 

Coefficients in the effective Hamiltonian, eq. (f|), for A% IS = 200 MeV and m b = 4.8 GeV, 
at a scale \i = m c = 1.5 GeV, in different renormalization schemes. See text for further 
details. 



10 



decay channel 


10 2 x 
(exp.) 


10 2 x B.R. 
(th.) 


10 3 x acp 
cos 5 > 


10 3 x acp 
cos 5 < 


D+ -> 


< .53 


.21 


— 


— 


D+ -> 7T + 77 


.66 ± .22 


.78 


-1.5 ±0.4 


-0.7 ±0.4 


L>+ -> tt+t/ 


< .8 


.10 


0.04 ±0.01 


0.02 ±0.01 




.73 ± .18 


1.11 


1.0 ±0.3 


0.5 ±0.3 




< .12 


.20 


-2.3 ±0.6 


-1.2 ±0.6 


D + -> p+TT 




.39 


2.9 ±0.8 


1.5 ±0.8 


£>+ -> p+77 


< 1.0 


.12 


— 


— 


D+ -> p+r/ 


< 1.4 


.07 


— 


— 


L> + -> + 


< .6 


.07 


— 


— 


£>+ -> 07T+ 


.60 ± .08 


.36 










1.85 


-0.9 ±0.3 


-0.5 ±0.3 




.47 ± .09 


.36 


2.8 ±0.8 


1.4 ±0.7 



TABLE 2 



Branching Ratios and CP-violating decay asymmetries for D + Cabibbo forbidden decays. 
Experimental data and 90% c.l. upper bounds taken from ref. |TTJ . 
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decay channel 



10 2 x B.R. 

(exp.) 



10 2 x B.R. 
(th.) 



10 s x acp 
cos 5 > 



10 3 x acp 
cos 5 < 



.12 



-0.5 ±0.2 



-0.2 ±0.1 



D+ -> 



.004 



2.3 ±0.7 



1.2 ±0.6 



.80 



1.0 ±0.3 



0.5 ±0.3 



if U 7T 



< 



.6(°) 



.48 



-1.5±0.5 



-0.8 ±0.4 



D+ -> 



1.59 



0.6 ±0.2 



0.3 ±0.2 



.15 



1.0 ±0.3 



0.5 ±0.3 



K+ 



.03 



-5.5 ±1.6 



-2.7 ± 1.4 



+ 



< 



. 2 W 



.18 



0.14 ±0.04 



0.07 ±0.04 



.03 



-5.8 ±1.7 



-2.9 ± 1.5 



.05 



-3.5 ±1.0 



-1.7 ±0.9 



.01 



3.1 ±2.3 



-4.1 ±2.1 



K*°7T + 



.23 



-2.6 ±0.8 



-1.3 ±0.7 



TABLE 3 

Branching Ratios and CP-violating decay asymmetries for Df Cabibbo forbidden decays 
Experimental 90% c.l. upper bounds are taken: (a) from ref. [HI and (b) from ref. |T7j . 
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